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Abstract 

A scaling consideration of the Kondo lattices is performed with account of 
singularities in the spin excitation spectral function. It is shown that a non- 
Fermi-liquid (NFL) behavior between two critical values of the bare s — f 
coupling constant occurs naturally for complicated magnetic structures with 
several magnon branches. This may explain the fact that a NFL behavior 
takes place often in the heavy- fermion systems with peculiar spin dynamics. 
Another kind of a NFL-like state (with different critical exponents) can oc- 
cur for simple antiferromagnets with account of magnon damping, and for 
paramagnets, especially with two-dimensional character of spin fluctuations. 
The mechanisms proposed lead to some predictions about behavior of specific 
heat, resistivity, magnetic susceptibility, and anisotropy parameter, which can 
be verified experimentally. 
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I. INTRODUCTION 



Recently, a great experimental material has been obtained for anomalous /-systems 
demonstrating so-called non- Fermi-liquid (NFL) behavior Manifestations of the 

NFL behavior are unusual temperature dependences of magnetic susceptibility (x(T) ~ 
T~^,( < 1), electronic specific heat (C(T)/T is proportional to T - ^ or — InT) and re- 
sistivity (p(T) ~ < 2), etc. Such a behavior is observed not only in alloys where 
disorder is present (U^Yi^Pds, UPts^Pd^, UCus-^Pd^, CeCue^Auz, U^Thi^Beis), but 
also in some stoichiometric compounds, e.g., CeyNi 3 M, CeCu2Si 2 ,CeNi 2 Ge2 |3). The latter 
situation is most interesting from the physical point of view ||. There are a number of 
theoretical mechanisms proposed to describe the NFL state: two-channel Kondo scattering 
0J10' "Griffiths singularities" in disordered magnets |IO||lTl , strong spin fluctuations near 
a quantum magnetic phase transition |l2| , |il| , [i"4|| etc. Most of modern treatments of the 
NFL problem have a semiphenomenological character. The only microscopic model where 
formation of the NFL state is proven - the one-impurity two-channel Kondo model - seems 
to be insufficient, since important role of intersite interactions is now a matter of common 
experience 0. 

In the present paper we start from the standard microscopic model of a periodical Kondo 
lattice. Main role in the physics of the Kondo lattices belongs to the interplay of the on-site 
Kondo screening and intersite exchange interactions. This interplay results in the mutual 
renormalization of the characteristic energy scales: the Kondo temperature Tk and spin- 
fluctuation frequency To. We shall demonstrate that during the renormalization process 
"soft" boson branches can be formed, presence of singularities of spin spectral function 
being of crucial importance. Scattering of electrons by such soft collective excitations just 
leads to the formation of the NFL state (cf. fT4|). 

In Sect. 2 the renormalization group (scaling) equations are presented. In Sect. 3 we 
consider the antiferromagnetic (AFM) state with account of the spin-wave damping and the 
paramagnetic state with simple spin-diffusive dynamics. It turns out that in these cases a 
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NFL-like behavior (in a restricted temperature interval) is possible, especially in the case of 
quasi- two-dimensional (2D) spin fluctuations. In Sect. 4 we consider the general problem of 
singularities of the scaling function. In Sect. 5 we show that the NFL behavior up to lowest 
temperatures can be naturally obtained provided that we take into account magnon-like 
excitations in the case of a complicated spin dynamics. The excitation picture required is 
characteristic for real /-systems where several excitation branches exist. In Sect. 6 we discuss 
various physical properties and possible relation to experimental data. 



II. THE SCALING EQUATIONS 

To describe a Kondo lattice, we use the s — f exchange model 

H = ^cL c ka -1^2 SiCr a pct a Cip + </qS_ q S q + H a (1) 

kcr ia/3 Q 

where is the band energy, Sj and S q are spin-density operators and their Fourier trans- 
forms, I is the s — f exchange parameter, J q are the intersite exchange parameters, a are 
the Pauli matrices, H a is the anisotropy Hamiltonian which results in occurrence of the gap 
u in the spin- wave spectrum. In Refs. 0,0], the interplay of the Kondo effect and inter- 



site interactions was investigated by the renormalization group approach. The latter starts 
from the second-order perturbation theory with the use of the equation-of-motion method 
(within the diagram technique in the pseudofermion representation for the spin operators, 
such an approximation corresponds to the one-loop scaling). The relevant variables are the 
effective (renormalized) parameter of s — f coupling g e f(C) = —2pI e f(C) (C — > — is a flow 
cutoff parameter, p is the bare density of electron states at the Fermi level), characteris- 
tic "exchange" spin-fluctuation energy UJ ex (C), gap in the spin- wave spectrum uq(C), and 
magnetic moment S e f(C). In the magnetically ordered phase, ZJ ex is the magnon frequency 
cj q , which is averaged over the wavevectors q =2k where k runs over the Fermi surface; in 
the case of dissipative spin dynamics (paramagnetic phase) uJ ex is determined by the second 
moment of the spin spectral density. 
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Here we write down the set of scaling equations with account of the magnon damping 
1(C): 

dg ef {C)/dC = A,d\nS ef (C)/dC = -A/2 (2) 
d In UJ ex (C)/dC = -aA/2,d\nu (C)/dC = -6A/2 (3) 
d\n*f{C)/dC = -cA/2 (4) 

with 

A = A (C ,^ (C ),, o(C )) = ^(^,^),^), 

a = 1 — a for the paramagnetic (PM) phase, a = 1 — a', b = 1 for the antiferromagnetic 
(AFM) phase, a = 2(1 — a"),b = 2 for the ferromagnetic (FM) phase; a, a', a" are some 
averages over the Fermi surface (see Ref. ||15|| ). For the staggered AFM ordering in the 2D 
and 3D cubic lattices where 

d d 

J q = 2 Ji cos q, t + AJ 2 cos g« cos qj 

i=l i>j 

with Ji and J 2 being the exchange integrals between nearest and next-nearest neighbors 
(| Ji| ^> | J2 1), one can derive by using expansion in small q 

a'^2(rf-l)^|(exp(^kR 2 )) tk=0 | 2 (5) 

where R2 runs over the next-nearest neighbors, the electron spectrum is referred to the 
Fermi energy Ep = 0. Similar equations can be obtained for the Coqblin-Schrieffer model 



15| and in the case of anisotropic s — f coupling |16 
Using (|2])-([|) we obtain the explicit expressions 



uJex(C) = u ex exp(-a[g ef (C) - g]/2) (6) 
000(C) =oo exp(-b[g ef (C)-g}/2) 
S ef (C) = Sexp(-[g ef (C)-g]/2) 
7(C) = kg 2 ef (Cpexp(-c[g ef (C) - g}/2) 
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In the last equation of (|]), we have taken into account that the magnon damping is pro- 
portional to g 2 , 7 = kg 2 ZJ (the factor k is determined by the bandstructure and magnetic 
ordering), and the scaling equations should contain only renormalized quantities. 



III. THE SCALING BEHAVIOR IN THE PRESENCE OF DAMPING 



As demonstrated in Ref. [flql , the singularities of the scaling function 77 can result in 
occurrence of a NFL behavior in a restricted region due to fixing of the argument of the 
function 77 at the singularity during the scaling process, so that oj{C) — |C|, 

g ef {0 -g~2{£- X)/a, X = ln(£>/57), (7) 

(the scale D is defined by g e f(—D) = g,£ = In \D/C\). This region becomes not too narrow 
only provided that the bare coupling constant g = — 21 p is very close to the critical value g c 
for the magnetic instability {\g — g c \/ 9c ~ 10~ 4 -hl0~ 6 ). Here we consider the scaling process 
with account of a not too small magnon damping (only very small damping was introduced 
in Ref. (|IIJ to provide existence of the magnetic-non-magnetic ground-state transition at 
9 = 9c)- 

The scaling function 77 in the FM and AFM phases for simple magnetic structures reads 
(for simplicity, magnetic anisotropy is neglected in this Section) 

iKEW|ci,7/|ci) = - K- k < + h^f/c 2 )' 1 ) (8) 

\ v ' I t k =t k ,=0 

For an isotropic three-dimensional (3D) ferromagnet integration in (§) for 7 = const and 
quadratic spin-wave spectrum yields 

. 1 + z 2 (1 + x) 2 + z 2 1 + x z 1 — X z 

n{x,z) = — In arctan arctan (9) 

4x (l-x) 2 + z 2 z 1+x z 1 - x 

where x = ZJ ex /\C\, z = j/\C\. Note that last two terms in (^) play a role similar to that of 
the "incoherent" contribution to the function rj, which was treated in Ref. ||15| . 



For a 3D antiferromagnet integration in (|8[) for the linear spin-wave spectrum gives 
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v (x, z) = ^Re [(1 + iz) ln(l +x + iz) + (l- iz) ln(l + x- iz)\ (10) 
where we take into account the intersubband damping only, 

7q = n2I 2 S(J q - J Q )p 2 A q+Q , (11) 



The damping (11) can be put nearly constant at not too large q (the threshold value deter- 
mined by the AFM gap can be neglected due to formal smallness in J). In the 2D case we 
obtain in the same approximation 



V(x,z)= -47- 4- o (13) 



z/ 3 (x, 2;) 
z/ 4 (x, 2) + 

z/ 2 (x, z) = hi - x 2 + ^(1 -x 2 ) 2 + 4z 2 ] 



which modifies somewhat the result of Ref. |15| . 

The function A q determines, in particular, the factor k in (|5|). For a parabolic electron 
spectrum we obtain 



(4vr)- 1 (l -x 2 )- 1 / 2 , d= 2 



2X 



where x = q/2kp, 9(x) is the step function, z is the electron concentration (with both spin 
projections). 

Main Kondo renormalization of the magnon damping comes from its proportionality to 
the factor of S. Spin fluctuations can give correction to this factor, as well as for the magnon 
frequency |15| . For simplicity, we restrict ourselves in numerical calculations to the AFM 



case with 57 q /7 q = du^/u^ = 5S/S, so that c = 1 in (£|), (|j). The corresponding scaling 



trajectories are shown in Fig.l. One can see that, unlike Ref. ||15|| , the "linear" behavior, 
although being somewhat smeared, is pronounced in a considerable region of £ for not too 
small \g — g c \, especially in the 2D case. In the 3D case the linear region (|7|) is followed by 
a quasi-linear behavior with 
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9ef 



{t)~A{Z-\)MC) cc \C\ aA /\S ef {C) cc \C\ A ^ 2 , 



(15) 



where A < 2/ a. To investigate the latter behavior in more details, it is instructive to consider 
also the case of a paramagnet with pure dissipative dynamics. In the case of spin-diffusion 
behavior we have (cf. Ref. [ |15|] ) 



PM 



rj 



-i 



, uj = AVki 



(16) 



where T> is the spin diffusion constant. As demonstrate numerical calculations (see Fig.2), 
for g < g c the one-impurity behavior l/g e f(£) = l/g — £ is changed at £ ~ A by the behavior 



fljj) with 



A ^ [g ef (\)} 2 ^(0),g ef (\) ~ 0/(1 - \g), *(0) = ^(1) ~ 0.5 



(17) 



In the 3D case where 7] PM (x) = arctanx/x the quasi-linear NFL-like behavior g e f{0 takes 
place in a rather narrow region. However, in the 2D case we obtain 



1 + (1 + x 



2\l/2 
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(18) 



2(1 + x 2 ) 

and the NFL-like region becomes more wide due to a more slow decrease of t] pm (x) at 
x — > oo. Note that such regions are not observed for g > g c . 



IV. SINGULARITIES OF THE SCALING FUNCTION 

To get further insight into the NFL-behavior problem, we perform an analysis of singu- 
larities of the scaling function rj, which is more general than in Ref. [ I5[ . In the absence of 
damping we can write down 

77 (uJ ex /\C\,uj /\C\) = K (l - ool/C 2 )' 1 (19) 

q 

The singularities of rj correspond to the Van Hove singularities in the magnon spectrum and 
to the boundary points q = and q = 2kp- 

For q — > 2k p the magnon spectrum has the Kohn anomaly 



(q — 2kp) In \q — 2kp\, d = 3 
S^otT _ (20) 

lq 2 -Ak 2 F 6(q-2k F ), d = 2 



Taking into account the dependences ([14]), which hold qualitatively in a general situation, 
we obtain for v = C 2 — uJ 2 —>■ 

{In In If I, d = 3 
(21) 
6{v )v- 1 ' 2 , d = 2 

Note that the singularity of the form lnt> obtained for the 3D case in Ref. |15j (see also 
previous Section), is a consequence of a simplified "Debye" model of the magnon spectrum. 
In fact, such a dependence corresponds qualitatively to an intermediate asymptotics at 
approaching the singularity. In any case, a small damping of spin excitations should be 
introduced to cut the singularity. In the calculations below we put the damping parameter 
5 = 1/100 (see Ref. [13]). In fact, pronounced extrema of T], but not singularities themselves 



turn out to be important for the scaling behavior discussed below. 

For q — > we have A q oc q" 1 . Near the points of minimum (maximum) in the magnon 
spectrum we have — oc ±q 2 , and for v = C 2 — u) 2 m ->0we obtain 

± In \v I, d = 3 

V (v) oc <{ (22) 

T e(Tv)\v\-^ 2 , d = 2 

so that 1] — > — oo near the band bottom and rj — > +oo near the band top. The Van Hove 



singularities in the magnon band at u = uj c for q ^ yield weaker singularities in r){v) (\v 
for d = 3 and a finite jump for d = 2) and will not be treated below. 
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V. THE SCALING PICTURE AND NFL BEHAVIOR IN MANY-SUBLATTICE 

MAGNETS 

Using the results of the previous Section we can propose a rather realistic and universal 
mechanism of the NFL behavior. Suppose that the spin excitation spectrum contains several 
branches which make additive contributions to the function rj. 
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As a simple example we can consider a two-sublattice ferrimagnet with the localized- 
system Hamiltonian 

Hf = ^(J q S_ q S q + J q S'_ q S' q + J q S_ q S q ), (23) 

q 

the s — / exchange interaction being taken into account only at one sublattice (spins without 
primes). Similar to ([19]) we obtain 

1= E A^fl-^V 1 (24) 



q.,- 1.2 M<1 \ (1 



where 



5 q = {[5(J q - Jo) + S'(J' q - J' ) -{S+ S')J ] 2 - ASS'Jlf'\ (25) 
Wqi,2 = B n T |5( J q - Jo) - S"(J q - J^) + (S - S')J \ 

are the acoustical and optical modes. 
The dependence 

*(£) = E^t^i/^)^, (^/D)^] (26) 

i 

is shown in Fig. 3 for the case of two excitation modes in a 3D antiferromagnet (the expres- 
sions for the function r]i in the case of one mode with inclusion of anisotropy are given in Ref. 

The only property of the function ^, which will be important below, is the occurrence 
of the second zero with decreasing \C\ (or increasing £). This property follows immediately 
from existence of the "positive" singularity in \1/ near the maximum frequency of the lower 
branch and of the "negative" singularity near the minimum frequency of the upper branch 
(e.g., for u>i(q = 2kp) < uj 2 (q = 0)). One can expect that this is a general property of many- 
sublattice magnets. The singularities can be also connected with the crystal-field excitations 

As demonstrate both numerical calculations and analytical treatment, in some interval of 
the bare coupling parameter, g c \ < g < g c2 , the argument of the function \l/ becomes fixed at 
the second zero, C = Co, during the scaling process which is described by Eq.(§). This can 
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be illustrated for the simple case where a = b for all the modes (e.g., for an antiferromagnet 
in the nearest-neighbor approximation we have a = b = 1). On substituting (||) into (0) we 
obtain 



d{l/g ef )/dt = - a[g ef - g}/2) (27) 

Then we derive 

g ef {Z - oo) ~ (2/a)(£ - £o) - l/[^(£o)£ 2 ] (28) 

where £o ~ \sx{D/W) is the second zero of the function Note that the first zero does 

not work since the unrestricted increase of g e f corresponds to a decrease of the argument of 
the function \1> in (|27D , so that \I/ — > +0. According to we obtain 

^(C),^) oc |C|,S e/ (C) oc ICI 1 ^ (29) 

Within the approach used, the behavior (p^),(p5|) takes place up to C = (£ = oo). Al- 
though the one-loop scaling equations themselves may become invalid with increasing g e f, 
the tendency to the formation of the "soft" magnon mode seems to be physically correct. 
The scaling picture for three possible cases is shown in Fig. 4. One can see that the interval 



[<7ci j Ate] where the NFL behavior occurs is not too small, unlike Ref. fig . 

In a more general case, where a ^ b and the exponents in (|3]) differ for different frequen- 
cies, the linear C-dependence takes place only for the total characteristic frequency ZJ (e.g., 
for an anisotropic antiferromagnet with one mode we have uj 2 = Uq+qJ 2 x ), and the behavior 
9ef(0 an d S e f(£) is more complicated. As follows from (|5|), for the AFM state with small 
next-nearest exchange interactions the case of small \a — 1| and b = 1 is realized. 

VI. DISCUSSION OF PHYSICAL PROPERTIES AND CONCLUSIONS 

Consider the temperature dependence of the magnetic susceptibility. In the spin-wave 
region we have for an AFM structure with the wavevector Q 
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X = lim«^|^ q )> w=0 = (Jo - Jq)- 1 oc S/uj (30) 

One can assume that the spin-wave description of the electron-magnon interaction is ad- 
equate not only in the AFM phase, but also for systems with a strong short-range AFM 
order (e.g., for 2D and frustrated 3D systems at finite temperatures). Anomalous /-systems 
demonstrate indeed pronounced quasi-2D spin fluctuations, see, e.g., Refs. |14] , |18| . Using 
the scaling arguments we can replace U — > 57(C), S — > S e f(C) with |C| ~ T, which yields 

x(T)QcT-f,C=(o-l)/a (31) 

According to (|5]), the non-universal exponent A is determined by details of magnetic structure 
and can be both positive and negative. For a qualitative discussion, we can still use Fig. 4 and 
treat the difference perturbation. The increase of x{T — > 0) (which is usually called 

NFL behavior) takes place for a > 1 and, as follows from (El), corresponds to an increase of 
magnetic anisotropy parameter with lowering T (see Ref. ||16|| ). Such a correlation may be 
verified experimentally. 

The temperature dependence of electronic specific heat can be estimated from the second- 
order perturbation theory, C e i(T)/T oc 1/Z(T) where Z(T) is the residue of the electron 



Green's function at the distance T from the Fermi level (cf. Ref. |17|). Then we have 



C el (T)/T ex gif(T)S e f(T)/u ex (T) oc X (T) In 2 T (32) 

The dependence C e i(T)/T cx x{T) has been recently obtained experimentally for a wide 
class of NFL systems ||11|| . One can expect that the accuracy of the experimental data is 
insufficient to pick out the factor of In 2 T. Generally, the temperature behavior of magnetic 
characteristics (S and ZJ), which depend exponentially on the coupling constant, is decisive 
for our NFL mechanisms. The transport relaxation rate determining the temperature de- 
pendence of the resistivity owing to scattering by spin fluctuations in AFM phase is given 
by 



\ = J/ 2 5 2 (Jo - J Q )p((v k+Q - v k ) 2 ) tk = £ A q (-?£) 

T V F q~Q \ °^q/ 



(33) 



q~Q 
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Then we obtain 

1 

T 



oc g 2 ef (T)S ef (T)/uJ ex (T)T 2 oc T 2 C el (T)/T oc T 2 ^ (34) 



Considering electron-electron scattering as the main scattering mechanism one can expect 
another temperature dependence, namely 

- oc [p/Z(T)f T 2 oc T 2 [C el (T)/T] 2 oc T 2 ~ 2? (35) 

T 

It would be interesting to compare experimental dependences on one hand, of the resistivity 
and, on the other hand, of the susceptibility and specific heat, to choose between dependences 



([34]) and (|35|) . Validity of the relation ( p4|) would be a verification of the scenario for NFL 
state formation proposed here. 

The temperature dependences of specific heat, magnetic susceptibility and resistivity in 
the case of a NFL-like behavior considered in Sect. 3 differ from those discussed above by 
the value of £ = (a — I) A/ 2 (for the second "linear" region). 

To conclude, phenomenon of the NFL behavior seems to have a complicated nature, so 
that it is hardly possible to propose an unified picture for all the cases. The mechanisms 
considered above are not based on disorder effects, but describe naturally the NFL state in 
ideal crystals. At the same time, the damping makes the quasi-NFL behavior considered in 
Sect. 3 more pronounced and in a sense plays the role of disorder. 

The damping is not important for the NFL behavior mechanism considered in Sect. 5. 
This NFL picture is "true" (i.e. holds up to lowest temperatures) within the lowest-order 
scaling approach; treatment of higher-order corrections to the scaling equations would pro- 
vide additional information. Unlike previous phenomenological works, existence of peculiar 
long-range critical fluctuations near the quantum phase transition is not needed for this 
mechanism, but local reconstruction of electronic states owing to the Kondo effect is essen- 
tial, the concrete form of spin spectral function being of crucial importance. More detailed 
investigations of the NFL behavior for complicated spectral functions, in particular with 
account of incoherent contributions, would be also of interest. 
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Figure captions 

Fig.l. The scaling trajectories g e f{0 in isotropic 2D antiferromagnets (solid lines, g = 
0.154 < g c ,g = 0.155 > g c ) and 3D antiferromagnets (dashed lines, g = 0.139 < g c , g = 
0.140 > g c ) with k = 0.5, a = c = 1, A = \n{D/uj) = 5. 

Fig. 2. The scaling trajectories g e f(0 in 2D (solid lines) and 3D (dashed lines) paramag- 
nets with a = 1/2, A = 5. The bare coupling parameters are g = 0.135 and g = 0.145, higher 
curves corresponding to larger g. The values of g c in 2D and 3D cases are 0.148 and 0.152. 

Fig. 3. The scaling functions defined by ( ^6|) in the case of a 3D antiferromagnet 
with two excitation modes The parameters are Z\ = 0.4, z 2 = 0.6, ln(D/uJ2) = 4, 0J2/0J1 = 3, 
w ,i/57e*,i = 0.2^0,2/^,2 = 0.6 (57? = uJ 2 ex>i + wj f ). 

Fig. 4. The scaling trajectories g e f{0 in a 3D antiferromagnet with the parameters of 
Fig. 3 and a« = hi = 1. The bare coupling parameters are g = 0.158 ~ g C 2 (upper solid line 
with the asymptotics), g = 0.154 ~ g cl (lower solid line), and g c i < g = 0.156 < g c2 (dashed 
line) . 
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